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Charles  K,  Sollitt 


Large  scale  laboratory  tests  were  conducted  to  determine  the  stability  of 
armor  rock  covering  an  existing  sewer  outfall  for  the  city  of  Goleta,  California. 
The  testing  consisted  of  two  phases:  phase  one  modeled  the  existing  condition  of 
the  outfall  and  phase  two  modeled  the  outfall  with  proposed  armor  rock  to  provide 
additional  stability.  Wave  tests  were  conducted  in  prototype  water  dq)ths  ranging 
from  15  to  45  feet  at  scale  ratios  ranging  from  1:4.52  to  1:15.5.  Prototype  wave 
conditions  included  both  random  and  monochromatic  waves  with  periods  ranging 
from  14  to  22  seconds  and  wave  heights  ranging  from  4  to  25  feet.  Wave  data 
were  taken  utilizing  seven  resistive  wave  gages,  two  acoustic  current  meters  and  a 
sonic  wave  profiler.  In  addition,  the  test  runs  were  video  recorded  from  two 
underwater  and  one  above  water  locations.  Test  conditions  for  each  test  run  are 
presented  in  tabular  form. 

Wave  conditions  were  analyzed  by  employing  Fourier  analysis  to  determine 
sine  and  cosine  amplitudes  of  each  frequency  component,  which  were  interpreted 
to  separate  the  incident  and  reflected  waves.  Results  are  summarized  in  tabular 
form  at  both  model  and  prototype  scale.  Significant  hydrodynamic  properties  are 


presented  graphically  in  non-dimensional  form  and  compared  to  theoretical  or 
empirical  models.  Similar  trends  were  observed  in  both  monochromatic  and 
random  wave  tests.  Breaking  wave  heights  were  found  to  be  within  75%  of  the 
theoretical  maximum  wave  height.  Reflection  coefficients  were  found  to  decrease 
and  transmission  coefficients  increase  with  increases  in  relative  water  depth. 
Non-dimensional  velocities  were  found  to  decrease  with  increasing  wave  steepness 
and  increasing  relative  depth.  Similar  hydrodynamic  trends  were  noted  in  the 
existing  and  proposed  armor  configurations.  Selective  removal  of  fines  from  the 
rock  armor  occurred  during  the  wave  tests.  The  proposed  armor  experienced  less 
removal  of  fines  than  the  existing  armor. 
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GOLETA  SEWER  OUTFALL  STABILITY  TESTS 


1.0  INTRODUCTION 


1.1  Background 

This  report  summarizes  wave  experiments  conducted  to  evaluate  the 
stability  of  armor  rock  covering  an  existing  marine  sewer  outfall.  The  subject  of 
this  study  is  the  Goleta  Sanitary  District  ocean  sewer  outfall,  located  in  the  city  of 
Goleta,  California,  approximately  nine  miles  west  of  Santa  Barbara.  Figure  1.1 
indicates  the  location  of  Goleta.  The  outfall  extends  5800  feet  from  shore  towards 
the  Santa  Barbara  Channel.  It  is  constructed  of  39  inch  diameter  welded  steel 
pipeline  and  has  a  tar  and  concrete  coating  which  produces  an  overall  outside 
diameter  of  44  inches.  Dispersion  of  effluait  is  through  a  280  foot  diffuser 
section  into  85  feet  of  water,  relative  to  mean  lower  low  water  (MLLW).  The 
primary  outfall  section  of  concern  extends  approximately  3600  feet  in  water  depths 
ranging  from  13  to  47  feet  with  an  average  bottom  slope  of  0.00945. 

The  Goleta  outfall  was  constructed  in  1965.  Existing  armoring  on  the 
outfall  consists  of  one  layer  of  Class  D  armor  rock,  with  a  median  weight  of 
approximately  0,5  lb,  which  was  placed  during  construction  and  one  layer  of  Class 
C  armor  rock,  with  a  median  weight  of  approximately  5  lb,  which  was  placed 
over  the  original  armoring  in  1981.  The  specific  gravity  of  the  existing  armor 
rock  is  2.60.  The  total  existing  armoring  reaches  an  elevation  of  8  to  12  inches 


Figure  1.  ’  Southern  California  coast  adjacent  to 
the  Santa  Barbara  Channel 
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below  the  top  of  the  outfall  and  extends  downward  on  both  sides  at  a  slope  of  1:3. 
However,  unsupported  lengths  of  pipe  have  been  observed  by  inspection  divers. 

Wave  forecasting  by  Brown  and  Caldwell  has  indicated  that  breaking  waves 
will  probably  occur  on  the  outfall  in  depths  shallower  than  47  feet,  with  incident 
wave  rays  approaching  at  a  mean  angle  of  22.5  degrees  from  the  longitudinal  axis 
of  the  outfaU.  The  proposed  method  of  stabilization  of  the  outfall  is  to  place 
layers  of  Class  A  and  B  armor  rock  over  the  existing  armor.  Figure  1.2  is  a 
cross-section  of  the  existing  outfall  and  proposed  armoring.  The  estimated 
specific  gravity  of  the  proposed  armor  rock  is  2.69.  Weight  distributions  are  as 
indicated  in  Tables  1.1  and  1.2. 

1.2  Stability  Analysis 

A  first  order  design  for  the  armor  rock  on  the  Goleta  outfall  may  be 
achieved  utilizing  the  semi-empirical  methods  presented  in  the  U.  S.  Army  Corps 
of  Engineers  Shore  Protection  Manual  (1984).  The  submerged  berm  formed  over 
an  armored  outfall  may  be  approximated  as  a  rubble  foundation  for  a  caisson 
structure  if  the  standing  wave  height  effect  is  included.  Refer  to  Figure  1.3  which 
is  reproduced  from  Figure  7-120  in  the  Shore  Protection  Manual.  The  caisson 
structure  will  cause  a  perfect  reflection  over  the  rubble  foundation,  effectively 
doubling  the  incident  wave  height  and  the  associated  destabilizing  wave  induced 
velocities  and  pressures.  Thus,  a  design  wave  height  of  H  for  rubble  foundation 
without  a  caisson  would  be  equivalent  to  a  design  wave  height  of  H/2  for  a  rubble 
foundation  with  a  caisson  structure. 


— ■? 
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2  LAYERS 
CLASS  A 


EXISTING  ARMOR 


Figure  1,2  Cross-section  of  suggested  wave  protection  for  Class  A  and  B  rock 
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Table  1.1  Prototype  size  distribution 
of  Class  A  rock 


Table  1.2  Size  distribution 
of  Class  B  rock 
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Figure  1.3  Stability  Number,  N„  for  rubble  foundation  and  toe  protection 
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Consider  a  breaking  design  wave  in  a  water  depth  of  IS  feet  at  the  outfall. 
This  will  yield  a  shallow  water  breaking  wave  height  of  approximately 

=  0.8J,  =  Uft 

The  equivalent  loading  condition  for  a  rubble  structure  with  a  caisson  would  be 
produced  by  a  wave  height  of  6  feet,  doubling  to  12  feet  by  a  perfect  reflection  off 
the  caisson. 

The  height  of  the  rubble  structure  is  6.7  feet,  as  shown  in  Figure  1.2. 

Thus,  the  depth  to  the  top  of  the  rubble  foundation  is 

d,  =  15. -6.7  =  8.3 /r 


and 


^  H 

d,  ■  15 


This  provides  a  stability  number  from  Figure  1.3  of  approximately 


For  a  specific  weight  of 


at/  *  20 


5,  =  2.69 


and  a  unit  weight  of  rock 


=  5^7^  =  2.69  (62.4)  =  168  Iblft^ 


the  stable  rock  weight  is 
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W= 


N]{Srif 


168(6)3 
20  (2.69-1)3 


=  376  lbs 


This  relatively  crude  calculation  verifies  the  80-98%  class  interval  rock  size 
for  the  Class  A  design  rock  in  Table  1.1. 


1.3  Scope 

This  report  describes  a  model  study  of  the  stability  of  existing  and 
proposed  armor  rock  on  the  Goleta  Sanitary  District  ocean  sewer  outfall.  Large 
scale  model  testing  was  conducted  at  the  Oregon  State  University  O.H.  Hinsdale 
Wave  Research  Laboratory  in  Corvallis  Oregon.  The  project  was  performed 
under  contract  with  Brown  and  Caldwell  Consultants  of  Irvine  California;  the 
model  construction  and  testing  took  place  between  February  10,  1992  and 
F^ruary  22,  1992.  A  weighted  PVC  pipe  was  constructed  and  armored  initially 
to  model  the  existing  condition  of  the  prototype  outfall,  shown  in  Figure  1.4.  It 
was  tested  in  prototype  depths  of  IS,  30,  and  45  feet  and  subjected  to  breaking 
waves  at  prototype  periods  of  14,  16,  19,  22  seconds.  The  model  was  also 
subjected  to  random  waves  modeled  as  a  JONSWAP  spectrum.  Next,  the 
structure  was  configured  to  model  the  proposed  Class  A  and  B  armoring,  as  in 
Figure  l.S,  and  subjected  to  a  similar  array  of  wave  conditions.  It  was 
understood  that  model  testing  of  the  proposed  Class  A  and  B  armoring  was 
dqpoident  upon  failure  of  the  first  model. 


#8  RE-BAR 


Figure  1 .5  Cross-section  of  model  for  proposed  armoring 


Wave  and  current  measurements  were  recorded  to  quantify  the  resulting 
wave  field  in  the  area  of  the  model.  Video  observations  were  recorded  for  each 
test  run  from  one  camera  above  water  and  two  cameras  below  water  to  monitor 
stability  of  the  armor  rock.  Wave  data  are  analyzed  using  Fourier  analysis 
combined  with  the  Goda  method  (Goda,  1985)  for  incident  and  reflected  wave 
separation.  A  summary  of  wave  and  current  conditions  is  presented  in  tabular 
form  at  both  model  and  prototype  scale.  Dimensionless  graphical  results  are 
presented  and  discussed. 
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2.0  MODEL  DESCRIPTION 


2.1  Facilities 

The  model  tests  were  conducted  at  the  O.  H.  Hinsdale  Wave  Research 
Laboratory,  Oregon  State  University.  The  major  features  of  the  Laboratory  include: 

1)  a  two  dimensional  wave  channel,  342  feet  long,  12  feet  wide  and  15  feet  deep  with 
a  hydraulically  driven,  hinged  flap  wave  board,  used  for  wave  induced  force  and 
response  measurements,  2)  a  three  dimensional  wave  basin,  87  feet  long,  60  feet  wide 
and  S  feet  deep  with  a  30  segment  directional  wave  generator,  electrically  powered 
through  ball  screw  drives,  used  for  wave  diffraction  studies,  3)  a  circular  wave  basin, 
50  feet  in  diameter  and  5  feet  deep  with  a  16  segment  spiral  wave  generator, 
electrically  powered  through  ball  screw  drives,  used  for  coastal  circulation  and 
sediment  transport  studies,  4)  a  VAX  serv^  3400  and  two  VAX  stations  3100  with 
optical  communication  links  for  wave  goieration  control  and  64  channels  of  digital 
data  acquisititm,  5)  a  50,000  ft^  oivironmental  oiclosure  for  the  facilities  including  an 
1,875  ft^  elevated  control  room. 

This  study  utilized  the  large  two  dimensional  wave  channel  to  maximize  the 
scale  of  the  model.  The  wave  board  is  hinged  at  the  bottom  of  the  channel  in  an 
overall  depth  of  18  feet.  The  wave  board  is  servo-hydraulically  drivra  with  direct, 
digital  ccmtrols.  A  150  horsqwwer  electrical  motor  powers  a  3000  psi,  76  gpm  oil 
hydraulic  pump,  which  is  the  prime  mover  for  an  dght  inch  diameter  actuator.  The 
actuator  ram  has  a  stroke  of  ±30  inches  and  is  located  10  feet  above  the  channel 
floOT.  The  back  sitte  of  the  wave  generator  is  dewatoed,  reducing  the  power 


consumption  by  one-half.  The  hydrostatic  head  is  overcome  by  applying  nitrogen  gas 
pressure  to  the  back  side  of  the  actuator,  creating  static  equilibrium  between  the  still 
water  pressure  and  a  gas  spring.  The  sides  of  the  waveboard  are  effectively  sealed  to 
the  sides  of  the  wave  channel  via  a  plastic  wiping  seal,  sliding  on  stainless  steel 
cladding  which  is  epoxied  to  the  concrete  channel  walls.  A  photograph  and  diagram 
of  the  wave  channel  appear  in  Figures  2. 1  and  2.2. 

Two  feedback  loops  are  used  to  control  the  waveboard.  The  primary  loop  is  a 
di^lacement  control  which  minimizes  the  error  between  the  measured  waveboard 
position  and  the  computed  position,  the  latter  based  upon  a  linear  wave  theory  solution 
for  the  wavd)oard  transfer  function.  A  secondary  loop  measures  the  wave  profile  at 
the  center  of  the  waveboard  and  corrects  the  board  velocity  to  yield  the  desired  wave 
profile.  This  loop  provides  the  capacity  for  active  reflected  wave  cancellation. 


Figure  2.1  O.H.  EQnsdale  Laboratory  2-dimensional  wave  diannel 


PLAN  VIEW 
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Figure  2.2  Plan  view  and  profile  of  the  2-dimensional  wave  channel 
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Both  monochromatic  and  random  waves  are  produced  with  this  wave  generator, 
over  a  useful  period  range  of  1  to  10  seconds.  Breaking  waves  up  to  5.0  feet  high  are 
generated  in  water  depths  of  11.5  feet,  for  wave  periods  of  4.0  seconds  or  less. 

Larger  waves  are  limited  in  height  by  the  stroke  of  the  wave  generator,  however,  long 
wave  breaking  can  be  achieved  by  shoaling  the  waves  with  a  movable  false  bottom. 

2.2  Scaling 

The  extreme  design  wave  condition  for  this  test  series  is  a  breaking  wave  at 
water  depths  of  45,  30  and  15  feet.  Given  the  desired  depth  and  wave  conditions  and 
the  physical  constraints  of  the  wavemaker,  it  was  determined  that,  with  shoaling,  a 
maximum  scale  of  1:10  was  possible  for  producing  a  breaking  wave  at  a  prototype 
depth  of  45  feet.  The  size  of  the  channel  allows  for  large  scale  model  testing  with  a 
minimum  distortion  in  viscous  characteristics  of  the  flow  field  relative  to  the  prototype 
condition.  As  long  as  the  armor  unit  Reynolds  number  of  the  model  exceeds  2x10^, 
the  error  in  viscous  effects  of  the  fluid  on  the  model  versus  the  prototype  is  less  than  3 
percait,  as  indicated  in  Figure  2.3  (SoUitt,  DeBok,  1976).  This  corresponds  to  a  one 
pound  armor  unit  and  a  four  foot  design  wave  height.  Smaller  armor  weights  and 
smaller  wave  heights  yield  smaller  Reynolds  numbers  with  a  corresponding  reduction 
in  the  apparent  model  stability.  This  provides  a  conservative  estimate  to  the  design 
condition  by  und^  predicting  the  no  damage  wave  height. 

Compressibility  of  the  liquid  is  of  no  consequence  in  this  application  and 
surface  tension  may  also  be  ignored  because  of  the  scale  of  the  waves.  Reynolds 
similarity  is  achieved  by  the  large  scale  of  the  model.  Because  gravity  is  the  restoring 
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Figure  2.3  Stability  correction  for  scale  effects  as 
a  function  of  Reynolds  Number 


\5 


force,  dynamic  similitude  may  be  characterized  by  maintaining  equality  of  the  Froude 
number  during  scaling.  The  Froude  number  maintains  the  ratio  of  inertial  forces 
(mVAV/AZ)  to  gravitational  forces  (mg)  or: 

F=(mVAVIAZ)lmg 


where,  m  =  mass  of  the  fluid  parcel 
g  =  gravitational  acceleration 
V  =  characteristic  velocity 
AV/AZ  =  velocity  gradient,  scaled  as  V/1 


thus. 


F=iVVt)fg^vVgl 


and. 


or, 


(w>/ir y^HviyfgT 


Givoi  that  the  Froude  number  remains  constant  and  a  length  scale  factor  (SF)  is 
characterized  by  the  equation: 


m  p 


thoi  the  following  scale  relationships  apply: 
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Length 

Time 

Velocity 

Weight 

2.3  Bottom  Profile 

Figure  2.4  is  a  profile  of  the  two-dimensional  wave  channel,  as  configured  for 
the  test.  Since  the  prototype  bottom  slope  over  the  3600-foot  section  of  outfall  under 
consideration  is  less  than  1:100,  the  model  test  section  was  placed  on  a  level  false 
bottom  surface.  The  elevated  felse  bottom  was  required  to  allow  shoaling  of  the 
waves  from  the  wavemaker  to  the  model.  Shoaling  is  necessary  to  provide  an  array  of 
desired  wave  characteristics  at  the  model,  given  the  limitations  of  the  wavemaker. 
Constraints  on  model  wave  height  include  wavemaker  stroke  and  velocity,  as  well  as 
wave  breaking. 

Target  wave  characteristics  at  the  wavemaker  are  calculated  by  determining  a 
desired  wave  at  the  model  and  then  shoaling  it  back  to  the  wavemaker  on  a  slope  of 
1:12  with  a  differential  height  of  10.75  feet.  The  level  false  bottom  section  was 
installed  using  12  foot  x  12  foot  x  8  inch  ndnforced  concrete  slab  sections  bolted  into 
place  and  caulked  between  sections.  It  extended  for  seven  sections  or  a  total  of  84  feet 
at  approximately  midway  in  the  channel.  A  beach  profile  on  a  slope  of  1:12  extended 
from  the  level  false  bottom  section  to  the  end  of  the  wave  channel,  providing  low 
reflection  for  a  wide  range  of  wave  periods. 
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2.4  Outfall  Pipe 

As  stated  previously,  a  maximum  scale  of  1:10  is  possible  for  producing  a 
breaking  wave  at  a  prototype  depth  of  45  feet.  However,  the  existing  armor  rock 
profiles  indicated  a  variation  in  gradation  with  depth.  Since  a  single  structure  was 
used  to  model  existing  conditions  at  each  of  the  three  indicated  depths,  the  scaling  was 
adjusted  for  each  depth  to  more  closely  represent  prototype  conditions.  A  scale  of 
around  1:12  allowed  flexibility  of  scaling  in  either  direction  keeping  within  the 
maximum  scale  limitation  of  1:10. 

A  prototype  44-inch  O.D.  outfall  was  to  be  modeled  at  a  scale  of  around  1:12 
which  gives  a  model  O.D.  of  3.67  inches.  In  consideration  of  the  commercial 
availability  of  materials,  a  3-inch  schedule  40  PVC  pipe  with  an  O.D.  of  3.5  inches 
was  chosen  for  the  model.  This  provided  a  model  to  prototype  scale  of  1:12.57.  The 
model  was  oriented  in  the  wave  channel  with  its  longitudinal  axis  at  an  angle  of  22.5 
degrees  fi’om  the  longitudinal  axis  of  the  channel,  giving  a  nominal  pipe  length  of 
28.97  feet.  Figures  2.5  is  a  diagram  of  the  pipe  section  orientation  and  Figure  2.6  is 
a  photograph  of  the  installed  pipe  section  prior  to  installation  of  the  armoring.  Since 
this  study  was  concerned  primarily  with  stability  of  the  armor  rock  and  not  the 
pipeline  itself,  the  model  pipe  was  weighted  with  six  12'foot  sections  of  No.  8  (1-inch) 
steel  reinforcing  bar  to  prevent  flotation  and  lateral  movement  under  wave  loading. 

2.5  Geologic  Materials 

The  existing  armoring  of  the  outfall  was  modeled  using  crushed  river  rock  with 


Figure  2.6  Installed  pipe  section  without  armor 


a  q)ecific  gravity  of  2.70.  The  rock  was  mixed  to  provide  a  model  distribution  which 
rq)resented  a  scaled  prototype  distribution.  A  sample  weight  distribution  of  the  model 
armor  rock  for  the  existing  condition  was  compared  with  a  prototype  sample 
distribution  taken  at  a  dq>th  of  30  feet  and  scaled  at  1:12.57.  Figure  2.7  is  a  plot  of 
the  model  and  scaled  prototype  weight  distributions  for  the  existing  condition.  As 
indicated  in  Figure  2.7,  the  model  rock  distribution  was  found  to  be  a  little  light  and 
the  model  scale  for  this  depth  was  adjusted  to  1:13.10,  which  represented  a  change  of 
2y)proximately  4.2  percoit.  The  same  model  was  used  at  a  scale  of  1:12.00  and 


7  Model  and  scaled  prototype  rock  weight  distributions 
for  the  existing  condition 
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1:15.50  tor  the  15-foot  and  45-foot  depths,  respectively,  based  on  variations  in 
prototype  rock  distributions  at  those  depths.  Scaling  the  wave  conditions  based  on  the 
armor  rock  scale  for  each  depth  produced  a  -5  to  -i-23  percent  error  in  modeling  the 
pipe.  It  was  considered  that  this  error  in  pipe  size  would  not  significantly  afiect  the 
stability  of  the  armor  rock. 

The  Class  A  armoring  was  modeled  using  crushed  basalt  with  a  specific  gravity 
of  2.77.  The  Class  A  rock  was  separated  by  sieving  with  a  Gilson  Model  TM-4 
shaker  and  sieve  sizes  of  1.5,  1.25,  1.0  and  .075  inches.  Mixing  by  weight  class 
intervals,  to  represent  a  scaled  prototype  distribution,  was  achieved  in  a  tumbling 
concrete  mixer.  The  process  is  illustrated  in  the  photograph  of  Figure  2.8.  The 
weight  distribution  of  the  mix  was  found  to  be  within  an  acceptable  range  of  variation 
of  the  prototype. 

Figure  2.9  is  a  plot  of  two  model  wdght  distributions  for  the  proposed  Class  A 
armoring  compared  to  target  minimum  and  maximum  limits.  The  Class  B  armoring 
was  modeled  using  pea  gravel  which  was  sieved  and  mixed  to  approximate  the 
prototype  Class  B  rock.  Figure  2.10  is  a  plot  of  the  model  size  distribution  compared 
to  target  limits  for  the  Class  B  rock  Although  the  Class  B  distribution  exceeded  the 
maximum  suggested  limit,  it  was  not  considered  that  this  would  significantly  affect  the 
outcome  of  the  armor  stability  tests.  The  proposed  Class  A  and  B  armor  rock  was 
scaled  at  12.57  with  no  variation  in  the  scale  with  dq}th. 

The  required  model  cross-section  of  the  existing  condition  was  constructed  by 
dropping  the  model  rock  through  an  18-inch  water  column  onto  the  pipe  section.  This 
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Figure  2.8  Rock  sieving  and  mixing  process 


procedure  is  illustrated  in  Figure  2.11.  The  model  armoring  was  not  compacted  or 
otherwise  disturbed  in  any  way  after  placement  except  that  the  toe  on  both  sides  of  the 
model  was  dressed  by  removing  excess  rock  that  fell  outside  the  limits  of  the  profile. 
No  high  spots  in  the  armoring  woe  adjusted  but  low  qx>ts  were  filled  in  with 
additional  rock,  providing  an  overall  armor  profile  that  tended  to  match  or  exceed  the 
required  armor  coverage.  Upon  completion  of  construction  of  the  first  model,  the 
wave  tank  was  drained  for  visual  inspection  prior  to  testing.  Figure  2. 12  is  a 
photograph  of  the  model  outfall  section  configured  to  the  existing  condition  of  the 


prototype. 


Figure  2.9  Model  weight  distributions  of  Class  A  armor  rock 
compared  to  target  minimum  and  maximum  values 


Figure  2.10  Model  size  distribution  of  Class  B  armor  rock 
compared  to  target  minimum  and  maximum  values 


Fig  2.11  Placement  of  model  rock  through  18-inch  water  column 


After  model  testing  of  the  existing  condition  was  complete,  the  model  was 
reconstructed  with  the  proposed  design  armoring.  The  required  cross-section  of  the 
second  model  was  constructed  by  dropping  first  the  Class  B  rock  and  then  the  Class  A 
rock  through  an  18-inch  water  column  onto  the  pipe  section.  The  wave  tank  was  then 
drained  and  a  plywood  template  was  used  to  manually  dress  the  sl(^  to  the  desired 
cross-section  of  Figure  1.2.  Figure  2.13  is  a  photograph  of  the  outfall  model  section 
configured  to  the  proposed  armoring. 


Fig  2. 12  Model  outfall  section  configured  to  the  existing  condition 


Fig  2. 13  Model  outM  section  ccmfigured  to  the  proposed  armoring 
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3.0  EXPERIMENTAL  PROCEDURES 


3.1  Overview 

Model  testing  was  conducted  over  a  two  week  period  from  February  10, 
1992  to  Fdjruary  21,  1992.  The  testing  was  divided  into  two  phases.  The  first 
phase  was  designed  to  model  the  existing  condition  of  the  Goleta  outfall  and 
determine  its  stability  under  various  monochromatic  and  random  wave  conditions. 
The  second  phase  was  designed  to  model  the  proposed  design  armoring. 

Prototype  sea  conditions  used  for  the  tests  were  prescribed  by  Brown  and 
Caldwell.  Wave  and  current  measurements  were  taken  on  the  seaward  and 
leeward  sides  of  the  model  to  quantify  the  actual  wave  conditions  of  each  test.  A 
binary  determination  of  stability  was  achieved  through  visual  observation  of  the 
model  during  each  test.  Video  cameras  recorded  each  test  from  one  above  water 
location  and  two  below  water  locations. 

3.2  Instrumentation 

Quantitative  data  recorded  during  each  test  included  wave  profile  and 
currrat  measurements.  The  wave  profile  was  measured  by  seven  resistive  type 
wave  gauges  spaced  along  the  east  side  of  the  wave  channel  adjacent  to  the  model 
and  one  sonic  wave  profiler  mounted  on  a  beam  at  the  center  of  the  channel  on  the 
seaward  side  of  the  mo(M.  The  resistive  wave  gauges  are  designed  and 
abstracted  at  the  Wave  Research  Laboratory  and  are  described  in  a  paper  by 
DiU>le  and  Sollitt  (1989).  Curroit  measurements  woe  taken  with  two  three- 
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directional,  Sensor  Data,  acoustic  current  meters  mounted  along  the  wall  of  the 
channel  on  either  side  of  the  model.  A  photograph  of  the  instrumentation  and 
underwater  video  camera  is  presented  in  Figure  3.1. 


Figure  3. 1  Acoustic  meter  and  video  camera  on  east  wall 


A  total  of  twelve  data  channels  were  used  to  input  wave  and  current 
information  into  the  laboratory’s  digital  data  acquisition  system,  one  for  each  wave 
profiler  and  two  (horizontal  and  vertical)  for  each  current  meter.  Data  was 
sampled  at  a  rate  of  16  data  points  per  second.  Table  3. 1  identifies  each 
instrumait,  its  corresponding  data  channel  and  its  location  in  the  wave  channel 
relative  to  the  center  of  the  model.  A  six  channel  pen  plotter  was  also  used  to 
visually  monitor  wave  and  curmit  measurements  during  tests. 
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Table  3.1  Instrument  Identification  and  Location 


Ch# 

Instrument 

Measurements 

Mounted/Side 

Ch# 

1 

Resistive 

Wave  Gage 

Wave  Profile 

-23 

ft . /East 

Ch# 

2 

Resistive 

Wave  Gage 

Wave  Profile 

-16 

ft. /East 

Ch# 

3 

Resistive 

Wave  Gage 

Wave  Profile 

-13. 

5  ft. /East 

Ch# 

4 

Resistive 

Wave  Gage 

Wave  Profile 

-l.C 

1  ft. /East 

Ch# 

5 

Resistive 

Wave  Gage 

Wave  Profile 

+25 

ft. /East 

Ch# 

6 

Resistive 

Wave  Gage 

Wave  Profile 

+32 

ft . /East 

Ch# 

7 

Resistive 

Wave  Gage 

Wave  Profile 

+35 

ft. /East 

Ch# 

8 

Acoustic  Current  Meter 

Ver.  Vel. 

-59 

in. /West 

Ch# 

9 

Acoustic  Current  Meter 

Hor.  Vel. 

-59 

in . /West 

Ch# 

10 

Acoustic  Current  Meter 

Ver.  Vel. 

+59 

in. /East 

Ch# 

11 

Acoustic  Current  Meter 

Hor.  Vel. 

+59 

in. /East 

Ch# 

12 

Sonic  Profiler 

Wave  Profile 

-60  ft. /Center 

(Minus  direction  is  toward  waveboard) 

Calibration  Date: 

Table  3.2  Wave  Gage 

Calibration  Constant: 
Feb.  10  Feb.  13 

Calibration 

feet  Tjer  volts 

Feb.  24  Average 

St.  Dv. 

Gage  No. 

1 

.4446 

.4508 

.4753 

.4569 

.0162 

2 

.4452 

.4377 

.4602 

.4477 

.0115 

3 

.4439 

.4393 

.4606 

.4479 

.0112 

4 

.4376 

.4320 

.4512 

.4403 

.0099 

5 

.4457 

.4265 

.4483 

.4402 

.0119 

6 

.4669 

.4462 

.4729 

.4620 

.0140 

7 

.4496 

.4306 

.4554 

.4452 

.0130 

Avg. 

.4476 

.4376 

.4606 

.4486 

.0115 

St.Dv. 

.0092 

.0087 

.0103 

.0082 
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In  addition  to  wave  and  current  records,  each  test  run  was  video  taped 
from  three  separate  locations  simultaneously.  Two  video  cameras  were  mounted 
below  the  waterline;  one  on  the  seaward  side  of  the  model  which  utilized  a  wide 
angle  lens  and  the  other  on  the  lee  side  of  the  model  which  utilized  a  telephoto 
lens.  The  third  video  camera  was  located  on  a  stationary  carriage  located  at  the 
center  of  the  channel  on  the  leeward  side  of  the  model  approximately  20  feet  away 
and  10  feet  above  the  still  water  level.  The  video  cameras  were  instrumental  in 
observing  the  stability  of  the  model  during  each  test  run.  Three  complete  and 
time  coordinated  video  tapes  of  each  test  run  were  recorded  including  audio 
comments  regarding  significant  characteristics  of  the  run.  The  video  tapes  are 
synchronized  with  each  other  and  the  digital  data  records  by  dubbing  WWV  time 
code  onto  audio  channel  2. 

Calibration  of  the  wave  gauge  instrumentation  was  performed  prior  to 
beginning  the  tests  to  ensure  accuracy  of  the  recorded  data  and  then  again  before 
the  second  phase  of  testing  and  upon  completion  of  all  tests.  Calibration  was 
achieved  by  filling  or  draining  the  channel  and  correlating  the  change  in  water 
depth  with  the  output  of  each  instrument.  Table  3.2  provides  calibration  constants 
for  each  resistance  wave  gauge  taken  during  each  calibration,  along  with  an 
avoage  calibration  constant  and  standard  deviation  for  each  gauge.  The 
calibration  constant  for  the  current  meters  is  10  volts  full  scale  per  meter  per 
second  of  velocity.  Figure  3.2  shows  the  location  of  the  instruments  used  in  this 
study. 


Figure  3.2  Diagram  of  instrumentation  in  the  wave  channel 
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3.3  Wave  Conditions 

The  majority  of  test  runs  were  conducted  with  monochromatic  or  simple 
periodic  waves.  The  repetition  provided  by  monochromatic  waves  allowed  for 
enhanced  judgement  of  the  stability  of  the  outfall  for  a  given  wave  condition.  A 
range  of  periods  was  chosen  based  upon  significant  wave  periods  that  were 
forecasted  at  the  prototype  location.  The  procedure  for  testing  the  outfall  stability 
in  monochromatic  waves  was  to  subject  the  model  to  waves  of  a  specific  period, 
increasing  the  wave  height  by  small  incremental  adjustments  to  the  stroke  of  the 
wavemaker  until  the  wave  was  breaking  at  or  near  the  model.  Monochromatic 
wave  tests  were  10  minutes  in  duration  and  the  period  of  the  waves  remained 
constant  during  the  run.  During  each  test  run,  the  stability  of  the  model  was 
observed  and  recorded  by  representatives  of  Brown  and  Caldwell. 

The  potential  for  damage  to  the  armoring  of  an  outfall  is  increased  in 
random  sea  conditions  where  a  short  series  of  very  large  waves  is  possible  during 
the  random  sequence.  Therefore,  random  waves  runs  were  included  as  part  of  the 
testing.  The  Jonswap  spectrum  was  chosen  to  model  the  random  sea  conditions. 
The  equation  for  this  spectrum  is; 

S(f)  =  oHl  Wltf)  eip[-1.25/(r,/)<]y’ 

where, 


33 


a  =  0.0624  /  [o.230  +0.03367-(0.185/(1.9  +  7))] 

0  =  exp[-(r^/-l)2  /  2<^] 
a  =  a^:f  <  fp,  a^:f  > 

^  0.07,  (Tj  t  0.09 

y  =  1  —  7  (mean  3.3) 

H,  =  significant  wave  height 
Tp  =  peak  wave  period 

The  significant  wave  height  and  peak  period  of  the  random  wave  sequence 
was  chosen  based  upon  forecast  prototype  design  wave  conditions.  The  procedure 
for  testing  the  model  in  random  waves  was  to  subject  it  to  a  random  wave 
sequence  varying  from  10  to  30  minutes  duration.  During  each  test  run,  the 
stability  of  the  model  was  observed  and  recorded. 

The  model  tests  w^  divided  into  two  phases,  designated  A  and  B.  Phase 
A  of  the  tests  modeled  the  existing  condition  of  the  outfall  for  breaking  and  non¬ 
breaking  monochromatic  waves  with  periods  of  14,  16,  19,  and  22  seconds  and 
with  random  waves  at  prototype  depths  of  IS,  30  and  4S  feet.  The  scale  of  the 
waves  varied  with  depth.  Phase  B  of  the  tests  modeled  the  proposed  additional 
armoring  on  the  outfall  for  breaking  and  non-breaking  waves  under  similar 
condititnis.  The  scale  during  phase  B  of  the  tests  remained  constant  for  all  test 
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runs  modeling  the  Goleta  outfall.  Figure  3.3  shows  monochromatic  waves 
brealdng  on  the  model  of  the  existing  condition  of  the  outfall.  Figure  3.4  shows 
the  model,  configured  with  the  proposed  armoring,  being  subjected  to  a  sequence 
of  random  waves. 

3.4  Summary  of  Test  Runs 

A  total  of  64  test  runs  and  three  calibration  runs  were  made  during  phase  A 
and  B  of  the  tests.  Each  run  was  identified  by  a  run  number  which  indicates  the 
phase  of  testing,  the  reciprocal  of  the  scale  factor  and  the  specific  run  designator. 
In  addition  to  tests  performed  specifically  for  the  Goleta  outfall  study,  tests  of  the 
San  Diego  Point  Loma  outfall  were  run  at  the  request  of  Brown  and  Caldwell. 

Test  runs  modeling  the  Goleta  outfall  are  identified  by  run  numbers  beginning 
with  A131,  A155,  A120  and  B126  where  the  numeric  code  represents  the 
reciprocal  scale  ratio  of  13.1,  15.5,  12.0  and  12.57,  respectively.  Test  runs 
which  modeled  prototype  conditions  which  were  non-specific  to  the  Goleta  study 
are  identified  by  scale  ratios  other  than  those  listed  above. 
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Figure  3.3  Monochromatic  waves  breaking  on  model  of  existing  condition 


Figure  3.4  Random  waves  breaking  on  model  with  proposed  armoring 
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4.0  RESULTS 


4.1  Overview 

The  quantitative  results  for  this  study  are  a  description  of  the  hydrodynamic 
environment  for  each  run  condition.  Wave  profiles  are  quantified  in  terms  of 
wave  period,  incident  and  reflected  wave  heights,  as  well  as  transmitted  and  beach 
reflected  wave  heights.  Equivalent  properties  for  both  random  and  monochromatic 
waves  are  reported  as  the  peak  wave  period  and  zero  moment  wave  height. 
Horizontal  velocity  amplitudes  of  the  currents  at  the  toe  of  the  slope  on  the 
seaward  and  leeward  sides  of  the  model  are  also  reported.  Results  are 
summarized  in  tabular  form  at  both  model  and  prototype  scale.  Graphical 
interpretations  of  dimensionless  results  are  also  provided. 

4.2  Analysis  Methods 

The  armored  outfall  pipe  creates  a  partial  blockage  of  the  water  column 
whidi  causes  a  fraction  of  the  incident  wave  energy  to  be  reflected  back  towards 
the  wave  generator.  Two  wave  profile  measurements  are  required  to  separate  the 
unknown  incident  and  reflected  wave  components.  This  is  accomplished  by 
spacing  two  wave  gauges  at  a  distance  between  5%  and  45%  of  the  wave  length 
and  performing  simultaneous  measuremoits.  Fourier  analysis  yields  the  sine  and 
coane  amplitudes  of  each  frequency  componrat  which  can  be  interpreted  to 
separate  the  inddoit  and  reflected  waves  (Goda,  1985).  This  method  is  applied  to 
bdh  monochromatic  and  random  waves. 
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The  zero  moment  wave  height,  is  the  height  of  a  single  wave 
component  which  has  twice  the  energy  of  the  measured  wave  system.  In  deep 
water  for  waves  of  low  steepness,  the  zero-moment  wave  height  is  equivalent  to 
the  significant  wave  height,  H,.  The  root-mean-square  wave  height,  H„bj,  has  an 
energy  level  that  is  equal  to  the  total  energy  in  a  random  sea  condition.  The  root- 
mean-square  wave  height  is  equal  to  the  zero  moment  wave  height  divided  by  the 
square  root  of  two. 


Thus,  a  random  sea  state  may  be  reported  as  an  equivalent  single  wave  in  terms  of 
the  peak  wave  period  and  either  the  zero  moment  wave  height  or  the  root-mean- 
square  wave  height. 

The  analysis  method  used  in  this  report  combines  Fourier  analysis  to 
separate  frequoicy  components  along  with  the  Goda  method  to  separate  incident 
and  reflected  waves.  The  total  energy  for  each  wave  system  is  obtained  by 
integrating  the  power  spectral  density  over  the  frequency  range  applicable  to  the 
measurement  scheme.  The  energy  of  the  incident  and  reflected  wave  trains  are 
interpreted  to  quantify  the  zero-moment  inddoit  and  reflected  wave  heights  on 
both  si(tes  of  die  structure.  The  ratio  of  the  reflected  to  incidoit  zoo  moment 
wave  height  provittes  an  estimate  of  the  effective  reflection  coeffidoit  for  the 
structure  and  beach.  The  peak  wave  periods  are  identified  at  the  peak  of  the 
measured  power  spectral  densities. 
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A  pair  of  wave  gauges  is  required'  to  compute  the  incident  and  reflected 
wave  characteristics  discussed  above.  With  three  gauges  each  on  the  seaward  and 
leeward  sides  of  the  model,  three  combinations  of  two  gauges  are  available  to 
separate  the  incident  and  reflected  waves.  This  provides  three  estimates  to  the 
v^ve  conditions  on  each  side  of  the  structure.  The  Appendix  to  this  report 
provides  all  three  estimates  of  incident  and  reflected  wave  conditions  for  each  test 
run  on  both  sides  of  the  structure.  Also,  an  analysis  of  Wave  Gauge  4  data  taken 
near  the  center  of  the  model  and  the  two  current  meters  is  provided.  Test  waves 
were  observed  to  continue  to  shoal  and  break  in  the  vicinity  of  the  model.  The 
observed  conditions  violate  the  stationary  and  spatially  independent  assumptions 
implicit  in  the  Goda  method.  Thus,  some  variation  in  the  results  between  different 
combinations  of  wave  gauges  is  expected.  Wave  induced  vdocities  of  the  toe  of 
the  seaward  and  leeward  slopes  of  the  structure  are  also  Fourier  analyzed.  No 
separation  of  incident  and  reflected  wave  components  is  possible  with  local 
vdodty  measurements.  Nevertheless,  the  equivalent  zero-moment  velocity  height, 
or  double  amplitude,  is  computed  for  each  record,  utilizing  a  method  analogous  to 
that  described  for  the  wave  profile  measuremoits.  The  zero  moment  vdocity 
double  amplitude  and  peak  period  for  the  two  current  measurements  are  also 
tabulated  in  the  Appendix. 

4.3  Tabular  Summary  of  Wave  Conditions 

Table  4.1  jmmdes  a  summary  of  the  modd  test  conditions  for  each 
separate  wave  test.  It  is  a  rdinence  log  for  the  analyzed  data  which  follows. 


Table  4.1  Data  Log  Summary 
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This  log  identifies  model  conditions  associated  with  each  numbered  run  including: 
scale  ratio,  date  and  time  of  experiment,  duration  of  data  sample,  whether  the  test 
utilized  random  or  monochromatic  waves,  the  water  depth  at  the  wave  board  and 
at  the  model,  the  input  wave  period  and  height,  the  oscillograph  observed  wave 
height  on  the  seaward  side  of  the  structure  and  the  fractional  equivalent  of  total 
wave  maker  span  required  to  achieve  the  model  wave  condition.  The  asterisk  at 
the  reported  random  wave  conditions  refer  to  the  input  peak  wave  period  and 
significant  wave  height. 

Table  4.2  provides  a  log  of  video  tape  numbers  associated  with  the 
recorded  underwater  and  above  water  video  documentation.  Note  that  each  video 
tape  is  spedtic  to  a  location  in  the  wave  channel  and  typically  includes  a  video 
record  of  several  consecutive  test  runs.  Therefore,  the  listed  test  run  number 
corresponding  to  a  particular  video  tape  represents  the  first  test  run  on  the  tape. 
The  locations  reported  in  the  table  are  to  be  interpreted  as  follows:  CALIB  refers 
to  calibration  observations  of  the  water  surface  rising  on  a  one-tenth  foot  division 
staff  on  the  side  of  the  wave  channel,  used  to  calibrate  the  wave  gauges;  Seaward 
refers  to  an  underwater  video  camera  attiu;hed  to  the  east  wall  of  the  wave  channel 
and  observing  the  seaward  slope  of  the  model  through  a  wide  angle  lens;  Leeward 
refers  to  an  underwater  video  camera  attached  to  the  west  wall  of  the  wave 
channel  and  observing  the  leeward  slope  of  the  model  through  a  telq)hoto  lens; 
Surfiu^e  refors  to  an  above  water  video  camera  placed  on  a  carriage  which  spanned 
the  width  of  the  wave  channel  and  observed  the  leeward  side  of  the  structure  and 
the  waves  interacting  with  the  structure. 
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Table  4.2  Video  Tape  Log 


Tape  No. 

Run  No. 

Tape  No. 

Run  No. 

1 

CAUB-1 

22 

B 126-01,  Leeward 

2 

CAUB-1 ,  cont'd 

23 

B1 26-01,  Surface 

3 

A1 31  -01 ,  Seaward 

24 

B 126-07,  Surface 

4- 

A131-01,  Leeward 

25 

B1 26-10,  Seaward 

5' 

A131-01,  Surface 

26 

B1 26-10,  Leeward 

6 

A131-08,  Surface 

27 

B126-15,  Surface 

7 

A131-11,  Seaward 

28 

B126-16,  Surface 

8 

A131-11,  Leeward 

29 

B126-16,  Seaward 

9 

A1 31 -13,  Seaward 

30 

B126-16,  Leeward 

10 

A131-13,  Leeward 

31 

B857-03,  Seaward 

11 

A131-13,  Surface 

32 

B857-03,  Leeward 

12 

A1 55-06,  Surface 

33 

B857-02,  Surface 

13 

A1 55-06,  Seaward 

34 

B750-01 ,  Surface 

14 

A 155-06,  Leeward 

35 

B750-01 ,  Seaward 

15 

A1 20-01,  Seaward 

36 

B750-01 ,  Leeward 

16 

A 120-01,  Leaward 

37 

Bl  13-03,  Surface 

17 

A1 20-01,  Surface 

38 

B1 1 3-04,  Leeward 

18 

CALlB-2 

39 

Bl  1 3-04,  Seaward 

19 

CALJB-2,  cont'd 

40 

CALJB-3 

20 

Armor  Rock  Const 

41 

CAUB-3,  cont'd 

21 

B1 26-01,  Seaward 

42 

Pump  Down  ^ 

The  Appendix  summarizes  the  analyzed  data  for  each  data  run  on  a  single 
page.  The  run  number  is  referred  to  as  "Test  Identification."  The  first  column 
summarizes  data  on  the  seaward  side  of  the  structure.  The  second  column 
summarizes  data  on  the  leeward  side  of  the  structure.  The  third  column 
summarizes  data  for  the  wave  gauge  at  the  center  of  the  structure  and  the  two 
current  meters.  At  the  top  of  the  first  two  columns,  the  model  water  depth  is 
listed  followed  by  the  "Data  Channels  us«I  to  compute  Coefficients"  and 
"Distance  between  channels  in  feet."  The  data  channel  numbers  refer  to  the 
numbers  in  column  one  of  Table  3.1.  The  distance  between  channels  refers  to 
wave  gauge  spacing,  which  is  the  difference  between  numbers  reported  in  colunm 
four  of  Table  3.1.  The  latter  is  required  input  to  Coda’s  analysis  for  separating 
the  incident  and  reflected  waves.  The  spacing  between  wave  gauge  pairs  used  in 
the  Coda  analysis  is  given  the  symbol  DL  and  is  reported  for  each  data  set.  Thus, 
in  column  one  of  the  Appendix,  combining  channels  1  and  2  corresponds  to  DL  = 
7.0  ft;  combining  channels  2  and  3  corresponds  to  DL  =  2.5  ft  and  combining 
channels  1  and  3  corresponds  to  DL  =  9.5  ft. 

For  each  of  three  wave  gauge  combinations,  a  summary  of  the  "Incident 
wave  energy"  and  "Reflected  wave  energy"  is  reported  in  the  Appendix.  In  each 
case,  the  raw  wave  energy  spectrum  is  smoothed  with  a  13  component  box  car 
averaging  scheme.  The  resulting  "Total  smoothed  energy"  is  listed  in  model  units 
of  ft^;  this  is  the  total  energy  integrated  under  the  power  spectrum  ci  r/e  and  is 
equal  to  the  variance  of  the  water  surface  profiles.  Energy  is  conserved  by  the 
smoothing  process.  The  "Maximum  smoothed  value"  is  the  peak  value  of  the 
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power  spectral  density  curve  in  model  units  of  fl?-sec.  The  "First  moment"  and 
"Second  moment"  of  the  power  spectral  density  are  related  to  shape  factors  of 
skewness  and  kurtosis,  indicating  asymmetry  and  peakedness,  respectively.  The 
zero  moment  wave  height,  "H,bo."  is  four  times  the  square  root  of  the  "T  otal 
smoothed  energy"  and  is  equivalent  to  the  significant  wave  height  for  deep  water 
waves.  It  is  the  wave  height  parameter  chosen  for  quantifying  the  water  surface 
profile  in  this  study.  The  wave  period  at  the  peak  in  the  spectral  density  curve  is 
"Tp."  In  some  instances,  very  asymmetric  monochromatic  waves  will  yield  at 
peak  at  a  sub-harmonic  or  super-harmonic  of  the  target  wave  period.  The 
"Reflection  coefficient"  is  simply  the  zero  moment  wave  height  of  the  reflected 
wave  divided  by  the  zero-moment  wave  height  of  the  incident  wave.  This 
information  is  repeated  in  the  Appendix  for  three  wave  gauge  spacings  on  the 
seaward  side  of  the  structure  (column  one)  and  three  wave  gauge  spacings  on  the 
leeward  side  of  the  structure  (column  2).  The  incident  wave  in  colunrn  two 
represents  the  wave  transmitted  beyond  the  structure.  The  reflected  wave  in 
column  2  represents  the  wave  reflected  from  the  1:12  sloped  concrete  beach. 

The  third  column  in  the  Appendix  presents  three  types  of  wave 
information.  The  first  group  quantifies  the  total  wave  energy,  incident  plus 
reflected,  at  Wave  Gauge  4  located  on  the  seaward  side  of  the  model,  on  the  east 
wall  opposite  the  model  centerline.  The  results  are  indicative  of  the  total  wave 
energy  at  the  structure.  The  quantities  presented  are  identical  in  character  to  that 
from  the  wave  gauge  pairs  except  the  reflected  wave  cannot  be  separated  from  the 
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incident  wave  with  a  single,  local  wave  measurement.  Thus,  no  reflection 
coefficient  is  quantified  from  Wave  Gauge  4, 

The  remaining  groups  of  data  in  column  3  of  the  Appendix  are 
interpretations  of  the  current  meter  measurements.  The  three  axis  sonic  current 
meters  were  located  six  inches  above  the  seabed,  approximately  one  foot  from  the 
toe  of  the  structure  on  seaward  and  leeward  sides.  Power  spectral  densities  of 
horizontal  and  vertical  velocity  components  were  resolved  from  these 
measurements.  Separation  of  incident  and  reflected  wave  effects  were  not  possible 
with  the  single,  local  velocity  measurements.  The  data  provided  is  analogous  to 
that  described  for  the  surface  wave  profile  measurements.  Note  that  the  zero- 
moment  velocity  height,  is  similar  to  the  zero  moment  wave  height.  It  is  a 
double  amplitude  representation  of  a  simple  harmonic  velocity  with  energy  equal 
to  twice  the  variance  of  the  actual  time  series.  It  could  be  interpreted  as  the 
double  amplitude  velocity  associated  with  the  significant  wave.  Identical  data  sets 
are  provided  for  both  the  horizontal  and  vertical  velocity  components  on  both  the 
seaward  and  leeward  sides  of  the  structure. 

Table  4.3  summarizes  significant  measured  hydrodynamic  parameters  from 
the  Appendix  at  model  scale  for  each  run  condition  identified  in  Table  4.1.  Only 
results  which  include  a  peak  period  within  25%  of  the  target  period  for  each  test 
run  are  summarized.  Where  multiple  estimates  of  incident  or  reflected  wave 
conditions  are  within  25%  of  the  target  wave  period,  an  average  of  the  tabulated 
values  is  computed.  If  the  peak  measured  period  is  not  within  25%  of  the  target 
period,  that  portion  of  the  summary  is  left  blank.  However,  the  full  data  summary 
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is  included  in  the  Appendix  for  all  run  conditions.  "Incident  wave"  characteristics 
are  summarized  from  the  seaward  side  of  the  structure,  column  1  in  the  Appendix. 
"Transmitted  wave"  characteristics  are  summarized  from  the  leeward  side  of  the 
structure,  column  2  in  the  Appendix.  "Model  Center"  present  results  from  Wave 
Gauge  4. 

Table  4.4  summarizes  the  measured  hydrodynamic  parameters  at  prototype 
scale  for  each  run  condition  identified  in  Table  4.1.  Values  reported  in  Table  4.4 
are  based  on  Table  4.3  values  which  are  scaled  by  Froude  number  as  discussed  in 
section  2.2  of  this  report.  Wave  height  and  water  depth  are  scaled  in  direct 
proportion  to  the  Iragth  scale  while  wave  period  and  velocity  are  scaled  in 
proportion  to  the  square  root  of  the  length  scale. 

4.4  Gr^hical  Summary  of  Wave  Conditions 

A  grEq)hic  presentation  of  significant  hydrodynamic  parameters  is  divided 
into  two  groups  corresponding  to  the  two  phases  of  testing.  The  two  phases  are 
the  existing  armor  condition  and  the  proposed  armor  condition.  For  each  phase, 
graphical  summaries  are  provided  of  the  wave  steepness,  reflection  and 
transmissitm  coefficients  and  dimensionless  wave  orbital  velocity.  The  results  are 
gr!q)hed  as  functions  of  the  dimoisionless  relative  dq>th,  i.e.,  the  ratio  of  the 
water  depth  to  the  linear  wave  theory  deep  water  wave  loigth.  The  latto* 
rqMesoits  a  rational  method  of  combining  water  dq>th  and  wave  pmod  as 
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17.420  17.086 


where  Lo  =  linear  wave  theory  deep  water  wave  length 
g  =  gravitational  constant 
T  =  wave  period 

Figures  4. 1  and  4.2  provide  graphs  of  the  non-dimensional  zero-moment 
wave  height  versus  non-dimensional  water  dq)th  for  various  depths.  A  theoretical 
breaking  wave  limit  (Dean  1974)  is  provided  for  comparison.  The  plotted  points 
determined  from  analysis  of  recorded  data  are  well  grouped  on  both  graphs  and 
follow  a  general  trend  similar  to  and  below  the  theoretical  limit.  This 
demonstrates  that  the  experimental  conditions  were  approaching  the  theoretical 
extreme  wave  environment.  It  should  be  noted  that  laboratory  wave  observations 
have  demonstrated  that  maximum  expoimental  wave  heights  are  approximately 
75%  of  the  maximum  theoretical  breaking  wave  height.  This  observation  is 
confirmed  by  Figures  4. 1  and  4.2.  The  figures  include  both  random  and 
monochromatic  test  results.  No  significant  trends  relative  to  absolute  values  of 
water  dq>th  are  jyjparent. 

Figures  4.3  and  4.4  are  graphs  of  the  non-dimensional  reflection  coefficient 
of  the  modd  versus  non-dimcmsional  depth  for  each  phase  of  testing.  Ahrens 
(1987)  provides  an  empirical  expression  for  wave  reflection  from  a  reef 


breakwater  as  follows: 


L 1  Wave  steepness  versus  dimensionless  depth  for  existing  condition 
(relative  depth,  h*,  is  non-dimensionalized  by  pipe  diameter) 


Figure  4.3  Reflection  coefficient  versus  dimensionless  depth  for  existing  condition 
(relative  depth,  h*,  is  non-dimensionalized  by  pipe  diameter) 


Reflection  coefficient  versus  dimensionless  depth  for  proposed  armoring 
(relative  depth,  h*,  is  non-dimensionalized  by  pipe  diameter) 
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where;  Ci  =  -6J1A  Cj  =  -0.0860 

Cj  =  -0.293  C4  =  -f-0.0833 

and;  d,  =  Water  depth  at  toe  of  breakwater 

Lj,  =  Airy  wave  length 
he  =  Crest  height  of  breakwater 
A,  =  Cross-sectional  area  of  the  breakwater 
F  =  Freeboard  of  the  structure 
H,no=  Incident  zero-moment  wave  height 
This  expression  is  plotted  as  Ahrens’  Equation  16  on  each  graph  for  comparison. 
In  each  phase  of  testing,  the  reflection  coefficients  obtained  from  analysis  of  the 
recorded  data  were  grouped  reasonably  well  but  were  considerably  higher  than 
those  calculated  from  Ahrens’  formula.  One  explanation  for  this  difference  is  that 
a  portion  of  the  transmitted  wave  energy  which  was  reflected  from  the  beach  slope 
at  the  end  of  the  wave  channel  was  again  transmitted  through  the  model  to  appear 
as  additional  reflected  wave  energy  from  the  model.  The  data  do  demonstrate  that 
longer  waves,  lower  h/Lo,  experience  greater  reflection.  This  is  consistent  with 
the  observed  behavior  of  other  coastal  structures  as  well. 

Wave  transmission  coefficients  for  the  model  are  similarly  presented  in 
Figures  4.5  and  4.6,  and  are  compared  to  Ahrens’  expression  for  wave 


Transmission  coefficient  versus  dimensionless  depth  for  existing  condition 
(relative  depth,  h*,  is  non-dimensionalized  by  pipe  diameter) 


Transmission  coefficient  versus  dimensionless  depth  for  proposed  armoring 
(relative  depth,  h*,  is  non-dimensionalized  by  pipe  diameter) 
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transmission  through  a  reef  breakwater.  Ahrens'  wave  transmission  coefficient 
(Ahrras,  1987)  is  plotted  as  Ahrens'  Equation  13  on  each  graph  and  is  given  by: 


where;  Cj  =  1.188  C3  =  0.529 

C2  =  0.261  C4  =  0.00551 

and;  h^  =  Crest  height  of  breakwater 

d,  =  Water  depth  at  toe  of  breakwater 
At  =  Cross-sectional  area  of  the  breakwater 
Lp  =  Airy  wave  length 
F  =  Freeboard  of  the  breakwater 
HnK)=  Incident  zero-moment  wave  height 
dso=  typical  dimension  of  the  median  stone 
Transmission  coefficients  were  calculated  from  the  results  of  analysis  using  the 
expression: 

K=H^(pransnatted)IH ^{incident) 

where  the  quantities  are  evaluated  in  Tables  4.3  and  4.4.  In  each  phase  of 
testing,  the  transmission  coefficients  obtained  from  analysis  of  the  recorded  data 
were  grouped  reasonably  well  and  were  in  a  range  to  be  expected  for  wave 
transmission  under  the  given  test  conditions.  In  both  cases,  Ahrens'  expression 
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plotted  near  unity  for  every  test  run,  which  indicates  that  his  expression  is  not  well 
conditioned  for  this  set  of  data.  The  recorded  data  imply  smaller  transmission  for 
long»  waves,  a  result  which  is  consistent  with  greater  reflection  being  observed 
for  longer  waves.  The  longer  waves  also  experienced  greater  losses  due  to 
breaking  on  the  structure,  thereby  reducing  the  transmission  coefficient. 

Horizontal  velocity  near  the  toe  of  the  seaward  side  of  the  model  is  non- 
dimensionalized  and  plotted  versus  non-dimensional  water  depth  in  Figures  4.7 
and  4.8.  Dean  (1984)  provides  theoretical  velocities  for  Case  A  (25%  of  the 
breaking  wave  height)  and  Case  D  (full  breaking  wave  height),  which  are  plotted 
for  comparison.  Since  the  velocities  determined  from  analysis  of  the  recorded 
wave  data  are  zero-moment  height  or  double-amplitude  velocities,  the  non- 
dimensionalized  values  from  Dean’s  tables  are  doubled  to  conform  to  the  same 
definition.  Similar  trrads  are  observed  between  the  theoretical  predictions  and  the 
measured  dimensionless  velocity.  In  general,  steq)er  waves  experience  greater 
decay  of  velocity  with  depth.  This  is  demonstrated  in  Case  D  lying  below  Case  A 
and  the  data  from  Gauge  4  (model  center)  occurring  near  the  minimum  of  the 
recorded  results.  Similarly,  the  transmitted  data  yields  the  larger  values  of 
dimensionless  velocity  for  minimum  values  of  wave  height. 

4.5  Selective  Removal  of  Rock 

During  the  frrst  phase  of  testing,  the  model  incurred  severe  wave  damage 
and  was  r^uilt  in  the  same  configuration  to  complete  that  phase  of  testing.  Post¬ 
test  rock  samples  were  takm  from  both  sides  of  the  model  and  from  the  rock  that 


Figure  4.7  Dimensionless  horizontal  velocity  versus  dimensionless  depth  for  existing  condition 


w. 

< 

Q 

0) 

0) 

C 

T3 

«> 

cS 

o 

o 

O 

c 

E 

(0 

(0 

Q> 

XJ 

a> 

•D 

CA 

C 

.  "c 
•  (0 

.  c 

1  1 

D  S 

■  Q 

t  t 

r-r 


Figure  4.8  Dimensionless  horizontal  velocity  versus  dimensionless  depth  for  proposed  armoring 
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was  washed  down  the  wave  channel  during  testing.  Figure  4.9  provides  post-test 
sample  rock  distributions  taken  upon  completion  of  the  first  phase  of  testing.  Note 
that  the  material  from  the  leeward  side  of  the  structure  is  considerably  heavier 
than  that  from  the  seaward  side,  suggesting  selective  removal  of  the  smaller  armor 
stone.  The  sample  of  material  washed-out  from  the  leeward  side  is  somewhat  less 
coarse  than  the  material  on  the  seaward  slope.  It  is  likely  that  the  fines  from  the 
washed-out  material  are  totally  removed  from  the  system  through  transport  up  the 
beach  and  settlement  through  the  gaps  in  the  beach  slabs. 

After  the  model  was  re-conf^gured  and  tested  during  the  second  phase  of 
the  project,  post-test  rock  samples  were  again  taken  from  both  sides  of  the  model 
and  from  the  rock  that  was  washed  down  the  wave  channel  during  testing.  Figure 
4.10  provides  post-test  sample  rock  distributions  taken  upon  completion  of  the 
second  phase  of  testing.  There  is  much  less  evidence  of  selective  removal  of  fines 
from  the  Class  A  armor  stone  than  from  the  existing  matnial.  The  leeward  slope 
material  and  the  washed-out  material  are  slightly  more  coarse  than  the  armor 
material  on  the  seaward  side. 


Figure  4.9  Existing  condition  post-test  rock  distributions 
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Figure  4, 10  Proposed  armoring  post-test  rock  distributions 
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5.0  SUMMARY  AND  CONCLUSIONS 


5.1  Test  Summary 

This  lepoit  summarizes  the  results  of  64  laboratory  tests  conducted  to 
determine  the  stability  of  armor  rock  covering  an  existing  sewer  outfall  for  the  city 
of  Goleta,  California.  The  tests  were  conducted  at  Oregon  State  University’s 
O.H.  Hinsdale  Wave  Research  Laboratory  in  Corvallis  Oregon.  The  testing 
consisted  of  two  phases:  phase  one  modeled  the  existing  condition  of  the  outfall 
and  phase  two  modeled  the  outfall  with  proposed  armor  rock  to  provide  additional 
stability.  Model  rock  was  obtained  from  local  quarries  to  reproduce  the  size 
distributions  of  existing,  and  proposed  armor  rock.  Wave  tests  w^  conducted  in 
prototype  water  depths  ranging  from  15  to  45  feet  at  scale  ratios  ranging  from 
4.52  to  15.5.  Prototype  wave  conditions  included  both  random  and 
monochromatic  waves  with  periods  ranging  from  14  to  22  seconds  and  wave 
heights  ranging  from  4  to  25  feet.  Wave  data  were  taken  utilizing  seven  resistive 
wave  gages,  two  acoustic  current  meters  and  a  sonic  wave  profiler.  In  addition, 
the  test  runs  were  video  recorded  from  two  underwater  and  one  above  water 
locations.  Model  test  conditions  for  each  wave  test  run  were  presented  in  tabular 
form.  Prototype  tabular  summaries  were  provided.  Non-dimensional  graphical 
summaries  were  presented  and  discussed. 
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5.2  Results  Summary 

Wave  conditions  were  analyzed  by  employing  Fourier  analysis  to  determine 
sine  and  cosine  amplitudes  of  each  frequency  component,  which  were  interpreted 
to  separate  the  incident  and  reflected  waves.  A  complete  set  of  the  results  of 
analysis  is  presented  in  the  Appendix.  Results  are  summarized  in  tabular  form  at 
both  model  and  prototype  scale.  Significant  hydrodynamic  properties  were 
presented  graphically  in  non-dimensional  form  and  compared  to  theoretical  or 
empirical  models.  Similar  trends  were  observed  in  both  monochromatic  and 
random  wave  tests.  Breaking  wave  heights  were  found  to  be  within  75  %  of  the 
theoretical  maximum  wave  height,  a  result  consistent  with  other  experiments. 
Reflection  coefficients  were  found  to  decrease  and  transmission  coefficients 
increase  with  increases  in  relative  water  dqjth.  Velocities,  non-dimensionaUzed 
with  respect  to  wave  height  and  period,  were  found  to  decrease  with  increasing 
wave  steepness  and  increasing  relative  depth.  Similar  hydrodynamic  trends  were 
noted  in  the  existing  and  proposed  armor  configurations.  Selective  removal  of 
fines  from  the  rock  armor  occurred  during  the  wave  tests.  The  proposed  armor 
experience  less  removal  of  fines  than  the  existing  armor. 
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